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Abstract 

Kostka numbers and Littlewood-Richardson coefficients play an essential role in the 
representation theory of the symmetric groups and the special linear groups. There 
has been a significant amount of interest in their computation (pQ, ^U|, [2], 0)- 
The issue of their computational complexity has been a question of folklore, but was 
asked explicitly by E. Rassart ^H]- We prove that the computation of either quantity is 
$R-complete. The reduction to computing Kostka numbers, is from the #P-complete 
problem [3] of counting the number of 2 x k contingency tables having given row and 
column sums. The main ingredient in this reduction is a correspondence discovered by 
D. E. Knuth j^j. The reduction to the problem of computing Littlewood-Richardson 
coefficients is from that of computing Kostka numbers. 

1 Introduction 

Let N = {1, 2, . . . } be the set of positive integers and Z> = N U {0}. Let A = (Ai, . . . , A s ) G 
N s , Ai > A 2 > ••• > X s > 1, H = (jii,...,Ht) G ^*> , v = (z/i,...,z/ u ) G Z| and 
a = . . . ,a v ) G N v , «i > • • • > a v > 1. The Kostka number and the Littlewood- 
Richardson coefficient c v Xa play an essential role in the representation theory of the symmetric 
groups and the special linear groups. Their combinatorial definitions can be found in Sec- 
tionEl The issue of their computational complexity has been a question of folklore. Recently, 
in POj, E. Rassart asked whether there exist fast (polynomial time) algorithms to compute 
Kostka numbers and Littlewood Richardson coefficients (Question 1, page 99). We prove 
that these two quantities are #P-complete (see Theorems C1I21); and thus answer his question 
in the negative under the hypothesis that a #P-complete quantity cannot be computed in 
polynomial time. 
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In Barvinok and Fomin show how the set of all non-zero K Xfl for a given \i can be 
produced in time that is polynomial in the total size of the input and output. They also give 
a probabilistic algorithm running in time, polynomial in the total size of input and output, 
that computes the set of all non-zero Littlewood-Richardson coefficients (f x given A and In 
j^j, methods for the explicit computation of the Kostka numbers and Littlewood-Richardson 
coefficients using vector partition functions are discussed. 

Combinatorially, the Kostka number K Xfl is the number of Young tableaux that have shape 
A and content \i ([5], page 25). The Littlewood-Richardson coefficient c Xa is the number of 
LR skew tableaux on the shape A* a with content v (this follows from Corollary 2, (v), page 
62 and Lemma 1, page 65 of [5]). Representation theoretically, K\p is the multiplicity of 
the weight \x in the representation V\ of the lie algebra sl r+ \(C) of the special linear group 
having highest weight A and c Xa is the multiplicity of V v in the tensor product V\ ®<c V a . 
They also appear in the representation theory of the symmetric groups (see chapter 7, [5]). 

While there are formulae for Kx^ and c Xa due to Kostant and Steinberg respectively (|3|, 
[2]), the number of terms is, in general, exponential in the bit-length of the input. These 
numbers have interesting properties such as, for fixed A, [i, a, u, K^xn^ and c^ XNa ([H]) are 
polynomials in N. Whether Kx^ > can be answered in polynomial time (see proposition |2J), 
and so can the question of whether c v Xa > 0, though the latter is a highly non-trivial fact 
estabilished by Ketan Mulmuley and Milind Sohoni jH] , and uses the proof of the Saturation 
Conjecture by Knutson and Tao [7j. This fact plays an important role in the approach to 
the P vs NP question [§] due to Ketan Mulmuley and Milind Sohoni. 

Let a = (ai,a 2 ) G Z> ,ax > a 2 and b = (&i, . . . bk) G Z> . We reduce the #P-complete 
problem ([!]) of finding the number |I(a, b)| of 2 x k contingency tables that have row sums 
a := (oi, a 2 ) and column sums b := (b 1 , . . . , bk), to that of finding some Kx^. We reduce 
this in turn to finding some c Xa , where A,//, a and v can be computed in time polynomial 
in the size of (a, b). The main tool used in the reduction to finding Kostka numbers is a 
correspondence found by Donald E. Knuth (jB]) between the set I(a, b) of contingency tables 
and pairs of tableaux having contents a and b respectively. 

2 Preliminaries and Notation 

A Young diagram ([5], page 1) is a collection of boxes, arranged in left justified rows, such 
that from top to bottom, the number of boxes in a row is monotonically (weakly) decreasing. 
The first two shapes in Figure 1 are Young diagrams. A filling is a numbering of the boxes of 
a Young diagram with positive integers, that are not necessarily distinct. A Young tableau 
or simply tableau is a filling such that the entries are 

1. weakly increasing from left to right across each row, and 

2. strictly increasing from top to bottom, down each column. 
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Figure 1: Left to right, the shapes A, a and the skew shape A * a. 

P and Q, in Figure 2, are Young tableau. A skew diagram is the diagram obtained from 
removing a smaller Young diagram out of a larger one. The third shape in Figure 1 is a skew 
shape. A skew tableau is a filling of the boxes of a skew diagram with positive integers, 
non- decreasing in rows, and strictly increasing in columns (see Figure 5). If the number of 
boxes in the i th row of a tableau, for 1 < % < s is A, and A := (Ax, . . . , A s ), it is said to have 
shape A. If the tableau houses fij copies of j for j < t and \x := (/ii, . . . , fi t ), it is said to 
have content /j,. Thus, in figure 2, P and Q have the same shape (5, 2), but contents (3, 2, 2) 
and (4, 3) respectively. 

We represent a Young diagram having s rows and Aj boxes in the i th row by A = (Ai, . . . , A s ). 
Given two shapes A and a, A * a is defined to be the skew-shape obtained by attaching 
the lower left corner of a to the upper right corner of A as in figure 1 (see page 60). 
size(A, /i) denotes the number of bits used in the description of this tuple of vectors. For 
A := (Ai, . . . , A s ), let |A| = J2i=i -V F° r vectors A, //, we say that A > /i if |A| = \fi\ and 
Vi, ^] Aj > YljKiVy I n addition, if A ^ //, we say A > \i. This ordering is called the 
dominance ordering. 

We call a tableau Littlewood-Richardson or LR, if, when its entries are read right to 
left, top to bottom, at any moment, the number of copies of i encountered is greater than 
or equal to the number of copies of i+ 1 encountered ([5], page 63). We denote the set of all 
(possibly skew) tableaux of shape A and content [i by T(A, //), and its subset consisting of 
all LR (possibly skew) tableaux by LRT(A, //). The Kostka number is the number of 
tableaux of shape A and content fi, i.e |T(A, //)| ((2J, page 25). The Littlewood-Richardson 
coefficient c u Xa is the number of LR skew tableaux of shape A * a of content u, i.e |LRT(A * 
a,v)\ (this follows from Corollary 2, (v), page 62 and Lemma 1, page 65 of |3]). 

3 The problems are in #P 

Valiant defined the counting class #P in his seminal paper [T2j, and proved that the com- 
putation of the permanent is #P-complete. The class #P is the class of functions / : 
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U ng N{0, l} n — > Z> , for which there exists a polynomial time turing machine M and a poly- 
nomial p such that (Vra G N), (Vac G {0, l}"),/(a;) = \{y\y G {0, 1} P ^ and M accepts 

The tableau shapes A, a and contents fx, v are described by vectors with integer coefficients. 
Therefore the number of boxes in a tableau might be exponential in the size of the input. 
Thus it needs to be estabilished that even though an object in T(A, fx) or LRT(A * a, v) may 
be exponentially large, it has a polynomial length description, and its membership (given 
this short description) can be verified in polynomial time. 

Proposition 1 The problem of determining Kx^ is in #P. 

Proof: 

Let A, fx be defined as in the Introduction. A tableau of shape A is fully specified given the 
number fiij of copies of i, that are present in the j th row, for all i < s and j < t. Thus its 
description is of length polynomial in size(A, fx), in fact it has length 0(size(A, fx) 2 ). Given 
the set {fJ>ij}i< s ,j<u this corresponds to the tableau obtained, by writing, in the i th row, from 
left to right, fXn copies of 1, fx i2 copies of 2 and so on for each i < s. To verify whether 
this filling has the content fx, we need to verify, that ^ (Xij = fXj for each j, which takes 
0(size(A, fx) 2 ) time. To verify whether this filling has the shape A, we need to verify, that 
fXij = Aj for each i, which again takes 0(size(A, fx) 2 ) time. To verify whether this filling is 
a tableau, we need to check that the entries in a column of this filling are strictly increasing 
from top to bottom. In other words, that, for each row i < s — 1 and each k < t — 1 

j<k j<k+l 

This can be done in time 0(size(A, fx) 2 ) if we maintained cumulative sums and calculated 

the next sum incrementally. Therefore the task of computing is in #P. 

□ 

Proposition 2 Given A and fx, whether or not > can be answered in polynomial time. 
Proof: 

Let A and fx be as in proposition [T] For any permutation a of the set {1, . . . , t}, let a(fx) be 
the vector (/V(i), • • • i/M*))- It is a known fact that = K X(T (^) (see jHj, page 26). Let a 
be a permutation such that Vz < t — 1, fx^ > fx a ( i+ i). For any fx, whose components are 
arranged in non-increasing order, it is known that K\p, > if and only if A > fx (see j^j, page 
26). Whether A > cr(fx) can be checked in time that is 0(size(A, fx)). Thus, whether or not 
K\fi > can be answered in time 0(size(A, fx) ln(size(A, fx)), which is the time it takes to 
find a permutation a that arranges the components of fx in non-increasing order. 
□ 

Proposition 3 The problem of computing c\ a is in #P. 
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Proof: 

Let A, a and v be denned as in the Introduction. Given S G LRT(A *a,u), we shall describe 
it as follows. Let the number of occurrences of j in the i th row of S be v^. S is described 
by (A, a) and {vij}i< s+v ,j< u - Conversely, given (A, a) and {vij}i< s+Vjj < u , we verify that it has 
shape A * a and content fx, by checking that 

1. Vi < v,J2j"ij = a h 

2. Vz > v + 1, J2j v ij = K-v, and 

This takes time 0(size(A, a, ^) 2 ). The description corresponds to the skew filling obtained, 
by writing in the shape A * a, left to right, in the i th row, copies of 1, v i2 copies of 2 and 
so on. Such a skew filling would be a skew tableau if and only if the entries of each column 
strictly increased from the top to the bottom. In other words, 

j<k j<k+l 

for each row i < s + v — 1 other than s and each k < u — 1. Thus, to verify from its 
description, that a given filling is an element of T(A * a, v) takes polynomial time. To verify 
that a specified filling is LR, we need to check, that, while scanning the entries of the tableau 
from left to right, top to bottom, at any instant, the number of copies of j encountered is 
> the number of copies of j + 1 encountered. In other words, that the following inequality 
holds for every row k < v + s and every j < u — 1 , 

i<k i<k+l 

Thus, to verify that a skew tableau is LR only takes time polynomial in the bit-length of its 

description. Thus, the problem of determining c\ a is in #P. 

□ 



4 Hardness Results 

Lemma 1 Given a = (ai, a 2 ) G Z> , a x > a 2 , and b = . . . , G Z> , let A = (|a|, a 2 ) 
and /x = (6i, . . . , b k ,a 2 ). Then, |I(a,b)| = K x ^. 

Proof: 
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Figure 2: An instance of the correspondence between I(a, b) and U^T(A,a) x T(A, b) for 
a= (4,3), b= (3,2,2). 

P Q p 

1 1 1 2 3 1 1 1 12 Discard Q 1112 3 

Construct Q 

2 3 2 2 2 3 



Figure 3: An instance of the correspondence between U^T(A, a) x T(A, b) and Uj vt>a T(A, b) 
for a = (4,3) and b = (3,2,2). 



The R-S-K (Robinson-Schensted-Knuth) correspondence (0, or [3] pages 40-41) gives a 
bijection between I(a, b) , the set of 2 x k contingency tables with row sums a and column sums 
b, and pairs of tableaux (T^T^) having a common shape but contents a and b respectively. 
In other words, we have a bijection between I(a, b) and U) V T(A, a) x T(A,b). A sample 
correspondence is shown in figure 2. 



Claim 1 For every shape X = (Ai,A2), such that that X j> a, there is exactly one tableau 
having shape A and content a. For any other shape A there is no tableau having shape A and 
content a. 



Any tableau with content a = (01,02) can have at most two rows, since the entries in a 
single column are all distinct. Further, 1 may be present only in the I s * row, because top to 
bottom, each column has strictly increasing entries. Therefore for a tableau to have content 
a, it is necessary that |A| = |a| and there be atleast a\ boxes in the I s * row. In other words, 

T(A, a) ^4> =^ A > a. 

Conversely, this condition together with a\ > 02, implies that a\ > A2. Therefore the filling 
in which the first a\ boxes of the top row contain 1 and all others contain 2 is a tableau (see 
Q in Figure 3). Since all the copies of 1 must be in the first row and must be in a contiguous 
stretch including the leftmost box, this is the only tableau in T(A, a). Hence the claim is 
proved. Thus there is a bijection between U^T(A, a) x T(A, b) and the set of tableaux of 
content b having some shape A > a. i.e, there is a bijection between U^T(A,a) x T(A, b) 
and U^ >a T(A,b). An example of this is provided in figure 3. Let us now consider the set 
U x>a T(X,b). 
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Figure 4: An instance of the correspondence between U^ >a T(A, b) and T(A,/i), where a 
(4, 3), 6 = (3, 2, 2), A = (7, 3) and p = (3, 2, 2, 3). 



Claim 2 Any tableau in U^ >a T(A, b) can be extended to a tableau of the shape A = (a± + 
a 2 , 02) by filling the boxes that are in A but not X, with k + 1. 27ms extension is a bisection 
between Uj it>a T(A, b) and T(A, fa). 



If there is a tableau of shape A and content a, Ai < |a| = ai + a 2 , and A2 < 02- A > a ==>- 
Ai > a 2 = A2. Therefore no two of the boxes in A which are not in A belong to the same 
column. Those of these boxes, that are present in a given row, occupy a contiguous stretch 
that includes the rightmost box. Therefore by filling them with k + 1 we get a tableau in 
T(A, /x). Conversely, given a tableau T in T(A,/i), deleting all boxes of T filled with k + 1 
gives a tableau in Uj vt>a T(A, b). These two maps are inverses of each other and hence provide 
a bijection between Uj lt>a T(A, b) and T(A, fa). Hence the claim is proved. 

An example of this correspondence has been illustrated in figure 4. Therefore, |I(a, b)| = 

I U A T(A,a) x T(A,b)| = | U A > a T(A,b)| = |T(A,/x)| = K x ,. 

□ 



Theorem 1 The problem of computing K\^, even when A has only 2 rows, is #P-complete. 



Proof: 

Computing is in by proposition ^ Now the result follows from Lemma Q because 

the computation of |I(a, b)| is known to be #P-complete (jl]). 

□ 



Lemma 2 Given A = (Ai,A 2 ) G Z> 0; Ai > A 2; and \i = (/^i,...,/^) G Z> , /et a = 
(ai,..., u>/iere (Vz)^ = J2j>i Vi> an d v = ("i, ■ ■ ■ ,ve), where Vi < £ — 1, u { = a { + fa, 

and is e = fa. Then K Xfl = c v Xa . 



Proof: 

c\ a is, by definition, |LRT(A * a, v) |, which is the number of LR tableaux on the skew shape 
A * a that have content v. The skew shape A * a consists of a copy of A and a copy of a, 
as in figures 1 and 5. For any skew tableau S of shape A * a, we shall denote by S\ a , the 
restriction of 5* to the copy of a and by S\\, the restriction of 5* to the copy of A. Thus, S\ a 
is a tableau of shape a and S\\ is a tableau of shape A. 
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Figure 5: An instance of the correspondence between T(A, fi) and LRT(A*a, is) for A = (7, 3) 
and fi = (3, 2, 2, 3), a = (7, 5, 3) and is = (10, 7, 5, 3). 



Claim 3 Let S G LRT(\ * a, is). For i < £ — 1, the i th row of S\ a must consist entirely of 
copies ofi. 



We proceed to show this by induction. Let i — 1. The rightmost entry of the I s * row of 
S\ a , is 1, by the LR condition. The entries of each row are non- decreasing, left to right. 
Therefore each entry of the 1 st row must be 1. Now assume for some r, 1 < r < £ — 1, that 
the claim is true for all i < r. If r = £ — 1, we are done, so let r < £ — 1. The entries of 
columns of S\ a are strictly increasing, top to bottom. By the induction hypothesis, the entry 
directly above a square in the (r + l) th row is r. Therefore, any element in the (r + l) th row 
must be > r + 1. But the rightmost entry of the (r + l) th row cannot be anything > r + 1, 
because this would violate the LR condition. The entries of each row are non- decreasing, 
left to right. Therefore each entry of the (r + l) th row must be r + 1. Hence, by induction, 
the claim is proved. 

Consequently, 5*1 a must have content v — a = /i. In other words, S'Ia G T(A, / u). Conversely, 
given any tableau T G T(A,/i), let S(T) be the skew tableau of shape A * a in which 
S(T)\x = T and the i th row of S^T)^ consists entirely of copies of i. While scanning S(T) 
right to left, top to bottom, the LR condition could not possibly be violated while on a 
square of S(T)\ a . By the time we begin scanning the rightmost box of the first row of 
S(T)\\, for any i < I — 1, the number of copies of i encountered is already /ii+i more than 
the number of copies of % + 1 encountered. Therefore the LR condition could not possibly 
be violated on any square of S(T)\\ either. Therefore S(T) G LRT(A * a, is). S(T)\\ = T, 
thus we have a bijection between LRT(A * a, is), the set of LR skew tableaux of shape A * a 
having content is and T(A, / u), the set of tableaux of shape A having content fi. Hence 
= |T(A,/i)| = |LRT(A * a, v)\ = c\ a as claimed. 

□ 



Theorem 2 The problem of computing c u Xa , even when A has only 2 rows is #P-complete. 
Proof: 

By Proposition El computing c v Xa is in #P. We have already proved in Theorem Q that the 



S 



computation of is #P complete. The result now follows from Lemma 121 
□ 



5 Conclusion 

We proved that the computation of Kostka numbers and Littlewood- Richardson coefficients 
is #P complete. The reduction to computing Kostka numbers was from the #P complete 
problem |4| of computing the number of contingency tables having given row and column 
sums. The problem of computing Kostka numbers was then reduced to that of computing 
Littlewood-Richardson coefficients. It is not known whether there exist Fully Polynomial 
Randomized Approximation Schemes (FPRAS) to compute any of these quantities. It would 
be of interest to find FPRAS to compute Kostka Numbers and Littlewood-Richardson coef- 
ficients. These are conjectured to exist by K. Mulmuley and M. Sohoni, jS] in their approach 
to the P vs NP question. 

6 Acknowledgements 

I wish express my gratitude to Ketan Mulmuley for suggesting the topic of this paper and 
for many valuable discussions. Many thanks are due to Ravi Kannan for informing me about 
[3]. I also sincerely thank Etienne Rassart, Rahul Santhanam, and Laszlo Babai for helpful 
comments. 



References 

[1] A. Barvinok and S.V. Fomin, Sparse interpolation of symmetric polynomials, Advances 
in Applied Mathematics , 18 (1997), 271-285, MR 98i:05164. 

[2] S. Billey, V. Guillemin, E. Rassart, A vector partition function for the multiplicities of 
sZ fc (C), Journal of Algebra, 278 (2004) no. 1, 251-293. 

[3] C. Cochet, Kostka Numbers and Littlewood-Richardson Coefficients, preprint (2003). 

[4] M. Dyer, R. Kannan and J. Mount, Sampling Contingency tables, Random Structures 
and Algorithms, (1979) 10 487-506. 

[5] W. Fulton, Young Tableaux, London Mathematical Society Student Texts 35 (1997). 

[6] D. E. Knuth, Permutations, matrices, and generalized Young tableaux, Pacific Journal 
of Mathematics, vol. 34, (1970), pp. 709-727. 



9 



[7] A. Knutson, T. Tao, The honeycomb model of tensor products I: Proof of the saturation 
conjecture, J. Amer. Math. Soc. 12 (1999), 1055-1090. 

[8] K. Mulmuley, personal communication, (2004). 

[9] K. Mulmuley and M. Sohoni, Geometric complexity theory, P vs. NP, and explicit ob- 
structions., Proceedings, International Conference on Algebra and Geometry, Hyder- 
abad, (2001). 

[10] E. Rassart, Geometric approaches to computing Kostka numbers and Littlewood- 
Richardson coefficients, Thesis for the Ph.D. degree in Mathematics, Massachusetts 
Institute of Technology (MIT), 2004. 

[11] E. Rassart, A polynomiality property for Littlewood- Richards on coefficients, Journal of 
Combinatorial Theory, Series A 107 (2004), no. 2, 161-179. 

[12] L.G. Valiant, The complexity of computing the permanent, Theoret. Comp. Sci., (1979), 
8:189-201. 



10 



